New insight in global QCD fits using Regge theory 
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In global QCD fits, one has to choose an initial parton distribution at Q 2 = Q 2 y I shall argue that 
the initial condition choses in usual standard sets is inconsistent with analytic 5-matrix theory. I 
shall show how one can combine these two approaches, leading to a Regge-compatible next-to-leading 
order global QCD fit. This allows one to extend the parametrisation in the low-Q 2 region. Finally, 
I shall discuss how it it possible to use the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equation to obtain information on Regge models at high Q 2 . 
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I. INTRODUCTION 



One of the most important results of QCD is the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution 
equation |j| which gives the evolution of the parton densities with the virtuality Q 2 of the photon. The kernel of 
these equation, i.e. the splitting functions, has been computed recently [Sly] at the three-loop level in perturbation 
theory (NNLO) . From the phenomenological point of view, the large number of experimental data coming from the 
highly-accurate measurement of the structure functions, e.g. at HERA 0,O|, allows for the determination of the 
{SJ parton densities. 

The technique is to choose a parametric distribution for the parton at an initial scale Q 2 ,, to obtain the Parton 
Distribution Functions (PDF) at all Q 2 using the DGLAP equation and to construct the physical quantities by 
convolution with the coefficient functions. One can then adjust the initial condition in order to reproduce the data 
as well as possible. This way of getting a standard PDF set, usually called a global QCD fit, has been applied many 
times with many updates by various teams |y, |g, III LLil • 

Beside this approach, it is well-known that one can reproduce the hadronic data using the analyticity properties 
of the S matrix. In this framework, the amplitudes are considered not as a function of the energy but in complex 
angular-momentum space. At high energy, the behaviour of the amplitude is then given by the leading singularities 
in the complex- j plane, corresponding to pomeron and reggeons exchanges. Therefore, in order to reproduce the 
data, one chooses a structure of singularity in the complex- j plane and adjusts their residues. Usually, we consider 
a pomeron contribution reproducing the growth of the cross-section at high energies, together with reggeon terms 
taking into account /- ,a- and w-meson trajectories exchanges. Within this framework, different models of pomeron 
are consistent with the present data As an example, one can consider the Donnachie-Landshoff two-pomeron 



model ^3 with one simple pole at j = 1.4 (F2 oc x 0A ) and a simple pole at j = 1.08 (F2 oc x ° 08 ). Throughout 
these proceedings, we shall concentrate on another choice known as the triple-pole pomeron model [Lj, where the 



pomeron singularity is a triple pole at j = 1, corresponding to a cross-section growing like log (l/x 



II. REGGE-COMPATIBLE GLOBAL FIT 



A. The initial condition problem 



The DGLAP equation, being a renormalisation group equation, gives you the the Q 2 -evolution of the PDF. This 
means that you still have to provide an initial condition at an initial scale Q\. Let us consider, e.g., the CTEQ 6M 
initial distributions. At small x, they obtain 

xq(x) oc £~ ' 30 and xg(x) oc x ' . 

This result suffers from two problems: first, the corresponding j-plane singularities are not seen in soft amplitudes. 
Secondly, quarks and gluons, while being coupled, do not have the same singularity structure. These two problems 
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TABLE I: Result of the global QCD fit together with the Regge small-Q 2 extension. 




FIG. 1: Valence quark, sea quark and gluon distributions for the NLO global QCD fit. 



are in contradiction with the fact that Regge theory requires jl4j all hadronic amplitudes to have the same j-plane 
singularities. 

B. Proposed solution 

Let us show how to solve this problem ^5|. We shall perform a global QCD fit starting at Qq = 5 GeV 2 with 
initial parton distributions compatible with Regge theory. In other words, one need to parametrise the initial parton 
distribution in agreement with S'-matrix theory. In order to obtain a complete PDF set, one need to parametrise 
the valence quarks u v and dy, the sea quarks u Sl d s , s s and c s , and the gluon density g. These expressions have 
to reproduce the pomeron and reggeon exchanges. Inserting a power of 1 — x in order to ensure that the parton 
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FIG. 2: Predictions of our fit for the longitudinal structure function Fl- 
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In these distributions, the log 2 (l/x), the log(l/x) and the constant terms correspond to a triple-pole pomeron exchange 
while the fourth one, with rj = 0.4 is a /-reggeon exchange. To build the initial conditions (|TJ, we have applied the 
following physical arguments 

• The pomeron is a flavour-singlet object. It therefore decouples from valence quarks and couples in the same 
way to all sea quarks. 

• The reggeon is expected to represent quark exchanges. We assume that it is not coupled to gluons. 

• The power of 1 — x defining the large- a; behaviour of the parton distributions is taken to be the same for all sea 
quarks. The exponent in the gluon distribution appears from a study of the DGLAP equation at large x |16|. 

• The mass effects has been introduced through the normalisation factors N s and N c . Due to the small strange 
mass, we expect N s » 1. N c is non-zero for initial scales Q 2 , larger than m 2 = 2 GeV 2 . Similarly, during the 
DGLAP evolution, the b quarks will be switched on at Q 2 — ml = 20.25 GeV 2 . 

• Quark number conservation fixes N u and Nd- The momentum sum rule is used to constrain C g . 

We have performed a global QCD fit with this initial condition at leading and next-to-leading order (in the MS 
scheme). We have included the data for the proton H. Isl Il7i IT3. Ilflj and deuteron 0, H3, EI structure functions, 
the ratio F% / F% [2lJ and the neutrino structure functions |22l |23| . The initial scale Qq has been set to 5 GeV 2 and 
we have imposed W 2 > 12.5 GeV 2 in order to cut the region where we expect higher-twist effects. We obtain, as 
presented in Table [I] a chi-square per data point slightly larger than 1. This is at least as good as the results reached 
by standard global fits. 

In Fig. 2]we show the valence and sea quark distributions. We also present the gluons obtained from our NLO 
fit compared with some distributions taken from various standard PDF sets at Q 2 — 5 GeV 2 . Our distributions are 
compatible with the standard PDF sets. 

Finally, one can consider predictions for the charm structure function F£ or for the longitudinal structure function 
Fl These are sensitive to the gluon distribution and, in addition, F| |25[ is a good test for the charm distribution. 
Our predictions are compatible with the data, as presented in Figs. [5] and [3 We also obtain a good description of 
the Tcvatron Drcll-Yan data with a if-factor between 1.3 and 1.4. 
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FIG. 4: Result of our fit for the NMC data for F 2 d (left) and F£/F% (right). In each plot, the high-Q 2 part is described by 
DGLAP evolution and the low-Q 2 data by Regge theory. 



III. SMALL-Q 2 DESCRIPTION 



One of the main advantages of restoring relevant analytic expressions for the initial distribution is that one can use 
26] Regge-theory techniques to describe the data for < Q 2 < Qq. We need to use the parametrisation and let 
the coefficients of the fit being Q 2 -dependent. For the Q 2 dependence, we have adopted for these dependences, the 
usual parametrisation 

where we have imposed matching with the parameters obtained from the global fit at Q 2 — Qq. With a few simplifi- 
cations, e.g., Q 2 D = Q% = Q% 3 = Q 2 Da , it gives a total of 13 parameters. 

At leading order, it is sufficient to consider extensions of the quark distributions. Due to the convolution with 
coefficient functions at NLO, one also needs small-Q 2 expressions for the gluon distribution and for the strong coupling 
constant a s . Since we expect our results to be relatively independent of this parametrisation, we have simply considered 
that the gluon density scales with Q 2 (g(x,Q 2 ) = Q 2 /Qog(x,QQ)) and that the running coupling stays constant for 
Q 2 < m 2 c . 

The x 2 resulting from this fit is presented in Table [I] together with the global QCD fit. We see that at NLO, one 
finally obtains a description of the DIS data, at all values of Q 2 , for W 2 > 12.5 GeV 2 with a chi-square of 1 per data 
point. 

To illustrate this, we have shown in Fig. 0]the final result of our fit compared with the and F% /F% measures 
from NMC. 



IV. REGGE THEORY AT HIGH-Q 2 

We have shown that it is possible to have a description of the DIS data based on Regge theory at small Q 2 and on 
DGLAP evolution at large Q 2 . However, it is well-known that DGLAP evolution generates an essential singularity at 
j = 1, which seems in disagreement with Regge theory. 



5 





FIG. 5: Sea quark (£) and gluon (xg) distributions at high Q 2 . The solid (resp. dashed, dotted) line corresponds to the fit 
starting at Qq — 10 GeV 2 (resp. 75, 550 GeV 2 ). For example, the solid curve has a triple-pole pomeron behaviour at Q 2 — 10 
GeV 2 and a essential singularity at other Q 2 values. 



It is proven that, in the small- a; region, additional resummation needs to be done if we want to recover a Regge-like 
high-energy behaviour. Our statement is that, even if their analytic structure are not physical, the DGLAP-obtained 
parton densities have to be treated as numerical approximations. Our aim is to show that, if we assume that DGLAP 
evolution approximates a triple-pole pomeron, we can extract the residues of the pomeron at large Q 2 from the 
DGLAP equation. In order to find the residues we may therefore adopt the following strategy |27j : 

1. choose an initial scale Qq at which one searches for the Regge residues, 

2. choose a value for the parameters in the initial distribution, 

3. compute the parton distributions for Qq < Q 2 < <5^ ax using forward DGLAP evolution and for Q 2 lin < Q 2 < Qq 
using backward DGLAP evolution, 

4. repeat and until the value of the parameters reproducing the Fi data for Q 2 > Q^ in and x < 0.15 is found. 

5. This gives the residues at the scale Q 2 and steps ^ to 0] are repeated in order to obtain the residues at all Q 2 
values. 

This technique allows us to obtain the residues of a Regge fit directly from QCD evolution, without having to 
postulate an analytic expression for the Q 2 dependence. It is also motivated by the fact that physics is expected to 
be independent on the choice of Qq. 

Applying this method between 10 and 1000 GeV 2 and for x < 0.15, we have obtained the <5 2 -dependent residues of 
the pomeron giving a fit to with x 2 /nop — 1.02. One may then check that the DGLAP result is a correct numerical 
approximation to a triple-pole pomeron model. In Fig. [31 we show the results of the DGLAP fit for 3 different choices 
of the initial scale. As seen from the upper part of the plot, in the case of the sea quark distribution, it is impossible to 
distinguish between the triple-pole pomeron behaviour (Q 2 — Qq) and the curves with a DGLAP-generated essential 
singularity (Q 2 ^ Qq)- This validates our initial argument. 
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The situation is not as clear for the gluon distributions (lower part of Fig. [SJ|. Since each of these fits reproduces 
equally well the F% data, the differences between them have to be considered as errors on the gluon distribution. We 
thus predict large errors on the gluon density at small Q 2 and small x, which may be of prime importance for LHC 
physics. 

V. CONCLUSIONS 

We have shown that it was possible to use Regge theory to constrain the parametrisation of the initial parton 
distributions in a global QCD fit. We obtain a standard PDF set at LO and NLO which reproduces the usual features 
of global QCD fits with the advantage of being consistent with analytic ^-matrix theory. 

This allows us to extend the fit in the small-Q 2 region giving a complete description of the data, for W 2 > 12.5 
GeV 2 , over the whole Q 2 range. 

Finally, we have discussed the compatibility of Regge theory and DGL AP evolution at large Q 2 . We show that one 
can use the DGLAP equation to extract the residues of the triple-pole pomeron and that this predicts large errors on 
the gluon distribution at small x and small Q 2 . 

In the future, it might be interesting to repeat this Regge-constrained global QCD fit with other Regge models and 
to combine the low- and high-Q 2 fits in order to allow for a determination of the initial scale for DGLAP evolution. 

The high-Q 2 study also invites one to look for QCD corrections to the DGLAP equation which stabilises the Regge 
behaviour. 
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